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ABSTRACT 

Dust formation in the winds of hot stars is inextricably linked to the classic eruptive 
state of luminous blue variables (LBVs) because it requires very high mass loss rates, M <; 
10~ 2 - 5 M Q /year, for grains to grow and for the non-dust optical depth of the wind to shield the 
dust formation region from the true stellar photosphere. Thus, dusty shells around hot stars 
trace the history of "great" eruptions, and the statistics of such shells in the Galaxy indicate that 
these eruptions are likely the dominant mass loss mechanism for evolved, Mzams ^ 40M Q stars. 
Dust formation at such high M also explains why very large grains {a mcLX <; 1 /im) are frequently 
found in these shells, since a max cx M. The statistics of these shells (numbers, ages, masses, and 
grain properties such as a max ) provide an archaeological record of this mass loss process. In 
particular, the velocities v s h e u » transient durations (where known) and ejected masses M s h e ii of 
the Galactic shells and the supernova "impostors" proposed as their extragalactic counterparts 
are very different. While much of the difference is a selection effect created by shell lifetimes 
fx (v shell V M shell)' 1 , more complete Galactic and extragalactic surveys are needed to demon- 
strate that the two phenomena share a common origin given that their observed properties are 
essentially disjoint. If even small fractions (1%) of SNe show interactions with such dense 
shells of ejecta, as is currently believed, then the driving mechanism of the eruptions must be 
associated with the very final phases of stellar evolution, suggestive of some underlying nuclear 
burning instability. 

Subject headings: stars: evolution; stars: massive; stars: mass-loss; supernovae: general; stars: 
winds, outflows; dust; extinction 



1. Introduction 

The winds of hot star s do not form dus t , as il lustra ted by the single uses of the word "dust" in the 



reviews of these winds by iKudritzki & Pulsl (120000 and IPuls et al.1 ((2008). The exceptions which prove 



the rule are the rare, dust forming WC stars, all of whi ch appear to b e binaries where dust forms due to 



the collision of the two stellar winds (see the review by ICrowtheri 120071) . and the relatively rare B[e] stars 



where dust is believed to form in a disk/dense equatorial wind surrounding the star (see the review by 
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Waters & Waelken s 1998). Dust formation is inhibited by the low particle densities and the harsh ultraviolet 
(UV) environment (e.g. Cherchneff & Tielenslll995l) . Cool star winds r eadily form du st, and there are 



extensive studies of dust formation in such environments (see the review by 



Willson 



2000). 



Yet it is clear the hot, massive stars can episodically form enormous quanti ties of dust and that this 
is related to the e ruptions of luminous blue variables (LBVs, see the reviews by Humphreys & Davidson 
1994 . 



Vink 



2009). LBVs are observed in three states: a quiescent, hot (T* > 15000 K) state, a cooler (T* 



7000 K) eruptive or S Doradus state of roughly the same bolometric luminosity but enhanced mass loss rates 
(M ~ 10 -4 to 10~ 5 M Q /year), and a similarly cool, "great" eruptive state of significantly higher bolometric 
luminosity and enormously enhanced mass loss rates (M ^> 10~ 2 -10~ 3 M Q /year). We will call these the hot 
(or quiescent), cool (or S Doradus), and (great) eruptive states, and we will refer to the ejected material from 
an eruption as a shell since the low duty cycles of eruptions produce relatively thin dusty shells of ejecta. In 
their (great) eruptive state, these stars can expel enormous amount s of material unde r conditions favorable 



to the growth of dust grains, as illustrated by the massive (~ 10M ffi , 



Smith et al 



2003 ), optical ly thick, dusty 



shell surrounding i] Carinae (see the reviews by Davidson & Humphreys 1997L Smith 2009 ). Indeed, it is 
likely that such phases represent the bulk of the mass loss from higher mass stars (M <; 40M(n) b ecause 



normal winds are inadequate to the task (Humphreys & Davidson! Il 984. 



recen t discovery of many 24/im shells surrounding hot stars by 



Smith & Owocki 



Wachter et al 



d20!0)and 



2006a). The 



Gvaramadze et al. 



(120101) further suggests that the phenomenon is more common than previously thought. 



Table 1 summarizes the p r opertie s of Galactic LBVs and LBV candidates mainly drawn from lHumphreys & Davidson 
(] 19941) and ISmith & Owockil (|2006ah . All the stars have substantial, M ~ 10" 5 M Q /year, relatively fast, 
Vqo ~ 200 km/s, present day winds that are not forming dust, as expected for hot stellar winds. However, 
all but one system is surrounded by a relatively massive, M^u ~ M & , slowly expanding, v s heii ~ 100 km/s, 
shell of dusty material, and in at least four cases, models of the shell appear to require surprisingly large 
maximum grain sizes, a max ^ 1 jum. We focused on these sources because most of these ancillary properties 
have been measured. Most of the mass estimates are based on assuming a dust-to-gas ratio = 0.01 and so 
could be underestimates. 

While Table 1 is certainly incomplete and subject to many selection effects, that it contains 13 objects 
means that shell ejections are an important or even dominant mass loss process for massive stars, as has 
been previously suggested in ord er to compensate for the steady downward revisions of t he mass loss rates 
in normal, hot stellar winds (e.g.. IHumphreys & Davidsonll 19841 ISmith & Qwockill2006al) . We can quantify 
this by estimating the number of dusty shells that should exist in the Galaxy given the rate of Galactic 
supernovae, r$N- For simplicity, we use a Salpeter initial mass function (IMF) and assume that supernovae 
arise from stars with initial masses M* in the range — SMq ^ ^ M up , where for now we will let 
M up — > oo. If eruptions occur in stars with M* > M erupt and there are an average of N erupt occurrences per 
star, then the eruption rate r erupt is of order 



'erupt 



0.1 



/4OM 



1.35 



NeruptrsN- 



(1) 



\M erupt 

The optical depth of a dusty shell of mass M s heii with visual opacity Ky ^ 100 cm 2 /g expanding at velocity 
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v shell is Ty = M s h e iinv /^v 2 hell t 2 so the shell will be detectable for of order 

( MshellKV \ 1/2 g . AA f°A ^ ( ' M shell \ ^ ( ™ ^ \ 

tsheii - — -J - 5700 — — — years. (2) 

\ 47rV shell T V J \ T V J V 10A W V v shell J 

Most shells should be seen near their maximum size, v s h e iit s heii — 0.4 pc, which is typical of the examples in 
Table 1. The total fraction of the stellar luminosity reradiated in the mid-IR is larger than ry because of the 
increased dust opacity in the UV. The expected number of Galactic shells is the product of the rate and the 
lifetime, 

As we see from Table 1, there are at least N s h e ii ~ 10 LBV stars surrounded by massive, dusty shells in the 
Galaxy, which means that the number of eruptions per star is 

N 9 (N sheU \(t v \ 1/2 f century 1 \ ( M erupt \ 135 ( 1OM \ 1/2 / v sheU \ 
Nm "-\\Q )W\) { r SN )\AmJ When) Ukm/sJ" W 

and the amount of ejected mass per star due to the eruptions is M tot = N erupt M s h e ii — 15M Q . A "normal" 
wind from a hot star with M norma i — 10 _5 M Q /year would have to operate continuously for over 10 6 years to 
equal the typical eruptive mass loss implied by the existence of even the well-studied Galactic shells. That 
Nempt > 1 is a lso consistent with the existen ce of multiple shells around some of the Galactic examples (e.g. 



G72.29+0.46, IJimenez-Esteban et al.ll2010f) . Note, however, that the overall duty cycle of the shell phase is 



low, since N erupt t s h e ii — 10 4 years as compared to post-main -sequence lifetimes of order 10 6 years. These 



estimates are broadly consistent with earlier estimates (e.g. Humphreys & Davidson! 1 1994, lLamerslll989l) 
but based on a different approach. 

Given such a large contribution to the mass loss history of massive stars, we need to understand the 
relationship between mass loss and dust formation around hot stars. We consider a parcel of fluid ejected 
in a wind of mass loss rate M and velocity v w ejected from a star of luminosity L* and temperature and 
examine the conditions under which dust can form in $2] Not surprisingly, the key variable is the mass loss 
rate. First, for stellar winds with velocities of order the escape velocities of massive hot stars, very high 
mass loss rates are needed for particle growth. Second, the dust formation region must be shielded from the 
hot stellar photosphere, which these high density winds can achieve by forming a pseudo-photosphere in the 
wind with a characteristic temperature of roughly 7000 K. Dust formation around hot blue stars is neces- 
sarily tied to very high mass loss rates, the classic LBV eruptive state and the formation of shells of ejecta. 
In $3] we discuss some implications of this model for dust formation, stellar evolution and supernovae. 



2. The Physics of Dust Formation in Stellar Transient Ejecta 



We consider the formation of dust grains of radius a comprised of Af atoms of average mass mo where 
47ra 3 Pbulk/^> = Nmo and py u n_ is the bulk density of the grain. We will use pbuik = 2.2 g/cm 3 (3.8 g/cm 3 ) and 
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Fig. 1. — Minimum mass loss rates M for silicate dust formation as a function of stellar temperature T*. 
The results for graphitic dusts are very similar. The lower left panel labels the regions. Dust can form 
in the region labeled "dust forms", above the photoionization ("ionized"), growth ("no growth") and pho- 
toevaporation ("photoevaporated") limits. The photoevaporation limits are for minimum photon energies 
of Eq = 7.5 (stronger) and 10 eV (weaker). The left, middle and right columns are for the three different 
assumptions about the apparent photospheric temperature. In the quiescent state (left) the photospheric tem- 
perature is the stellar temperature, T p hot = T*, in the S Doradus state (middle) the photospheric temperature 
of T p h () t = rnin (T* , 7000 K) is assumed to be determined by an expansion of the stellar photosphere that is 
uncorrected with the wind, and in the eruptive state (right) the apparent temperature is determined by the 
non-dust optical depth of the wind. The dotted contour in this column is the contour where T p h ot = 7000 K. 
The upper, middle and lower panels show the changes for stellar luminosities of L* = Lp^ot — 10 7 , 10 6 and 
10 5 L Q , respectively. The filled points are the present day properties of the systems from Table 1. Other 
parameters are set to M* = 2OM , v c = 1 km/s, X g = 0.005, T d = 1500 K, N = 7, p = 3.8 and m = 20m p . 
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mQ = \2m p (20m p ) for graphitic (silicate) grains. The smallest possible grain is the point where the inter- 
particle bond strengths shift from strong molecular bonds to weaker intra-molecular interactions, and we will 
generally require N >7 based on the number o f atoms in Mg?Si04. There are many prior tre atments of dust 



formation in (generally cool) stellar winds (e.g. Salpeterlll977 . Drainelll979 



Dwek 



1988 



Gail et al 



Kozasa et al 



1984), novae (see the 



1991) which contain 



review bv behrd[l988h and supernovae (e.g. Clayton! 1 19791 . 1 
most o f the b asic physical picture we use here. For some standard results in dust physics we will refer to 
Draind 020 1 lh as a reference source. For dust to grow, the medium must be largely neutral, sufficiently cool 



for growth to occur, have a high enough density for there to be an appreciable particle collision rate, and 
the grains must grow faster than they can be photo-evaporated by ultraviolet (UV) photons. We will assume 
that the first stage of particle formation, nucleation to form the smallest grains, simply occurs once the 
temperature is sufficiently low, and consider only the subsequent collisional growth of the grains. Making 
nucleation an additional bottleneck to dust formation will only strengthen our conclusions. 

The dust forms in a (time varying) wind which we can characterize by the mass loss rate M, the wind 
velocity v w , and the mass fraction of condensible species X g . Not all the condensible mass need condense 
onto grains, so the ultimate mass fraction of dust Xj < X g . The wind is produced by a star of luminosity 
L*, photospheric radius and effective temperature where L* = A-kR\oT^. At the time the dust is being 
formed, the star has luminosity L p j wt , radius R p i w t and temperature T p i rot , with L p \ wt = 4irR^ hot o-Tp hot . We 
consider dust formation in three physical states. First, we have the hot star in quiescence with L p \ wt = L* 
and T p hot = T*. Second, we have the S Doradus state, where L p h ot = L* and T p i wt = 7000 K^T*. Finally, we 
have the eruptive state where both L p h ot ^ L* and T p hot ¥ T* and we will use the optical depth of the wind to 
estimate Tphot- m many cases we can assume that the spectral energy distribution is simply a black body, but 
there are several areas where the differences between black bodies and stellar photospheres are important 
because line opaci t ies s uppress the ultraviolet emission. Where this is important, we will use the models of 
Castelh & Kuruczl J2003h . 



For the physics of dust formation, the most relevant velocity is that near the dust formation radius since 
it sets the particle density that determines the growth of the grains. For hot stars, dust formation occurs 
sufficiently far from the star that the wind acceleration should be largely complete and we can view the 
velocity as constant. We will scale the wind velocity by the surface escape velocity of the quiescent stars 



R* J \20MqJ V 







1/4 



10 4 K 



km/s, 



(5) 



as this is the typical asy mptotic velocity scale of radiatively accelerated winds (see the reviews by lKudritzki & Puis 



2000, 



Puis et al 



2008). Here we have scaled the stellar mass to M* = 20M Q under the assumption that the 



stars have undergone significant mass loss by the time of the eruption. Since v e oc T*, this introduces a strong 
stellar temperature dependence to dust formation and growth. In most of our results, the wind velocity is 
always scaled by the escape velocity of the quiescent star! Eqn. [5] yields an appropriate velocity scale for 
rj Carinae but may be somewhat high for many of the sources in Table 1 . For the more massive examples in 
Table 1 the low shell expansion velocities are almost certainly intrinsic because only a very high density in- 
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terstellar medium (ISM) can significantly slow the expansion of these massive shells of material^ Moreover, 
some of the systems have multiple shells where th e inner shells have lo w velocities b ut cannot be interacting 
with t he ISM or an older slower wind (e.g. P Cyg, Me aburn et al. Il 19961 : G79.29+0.46 ljimenez-Esteban et all 
20101) . Thus, it may be that some systems have asymptotic velocities significantly lower than those from the 
surface of the quiescent hot star - this can be approximated by replacing L* and I 1 * by L p ) rot and T p f, ot , which 
will change the scaling of the wind velocity to use the escape velocity from the photosphere of the transient. 

Luminous hot stars produce large numbers of ionizing photons, and dust cannot form in such a hot, 
ionized medium. For a pure hydrogen wind, the wind can reco mbine if the rate of production of ionizing 



photons is less than <2o = M aB/4Trv^m p R* (e.g. lFranssonlll982T) . The production of ionizing photons is 



Qo = T ^ r F(G = x- l ,E 1 /kT*,oo) 
kT* 



(6) 



where E\ = 13.6 eV, the dimensionless function is 



F(G, x ,xi) = 



P2 GF v dx 



15 



^F v dx 7T 4 J XQ exp(x)-l' 



Gx dx 



(V) 



where x = hv/kT and the limit is that of a black body. Thus, the minimum mass loss rate for the wind to 
recombine is 



M Z 



STrL*GM*m 2 p F 



1/2 



kT^as I 

3 1/2 /3 x 10" 13 cm 3 /s U 10 4 K Af* 



2.1 x lO" 3 ^ 1 / 



OtB 



10 6 L o J* 20Mcr 



1/2 



M /year. (8) 



For a black body, F ~ (15/7r 4 )(2+x(2+x)exp(-x)) where x = 158000/7*, making it a small number unless 
the star is very hot (F 1 / 2 = 0.0024 for a = 10 4 K black body), and the line blanketing of stellar atmospheres 
red uces it still further. Fig. Q] shows this photoionization limit ("ionized") on M for forming dust based on 



the lCastelli & Kuruczl (12003b model atmospheres - unless the star is very hot, ionization cannot prevent dust 
formation in dense winds. Even then, the photoionization limit is only important for the quiescent star. 

If the gas is relatively neutral, and so can carry out chemical reactions to form molecules, the tem- 
perature must be low enough to aggregate the molecules into grains. In general, the dust temperature is 
controlled by the radiation field because collisional time scales are very much longer (see below). We can 
divide the effects of radiation into the equilibrium temperature and stochastic heating of small grains by 



'Given the shell masses it is very hard to slow them down by large factors. Even slowing a IMq shell by a factor of two (from 
140 to 70 km/s) within an expansion radius of 0.5 pc requires an ISM density of order 10 2 cm 4 . A dense (M ^ lO~ 5 M0/year), 
pre-existing slow wind from a red supergiant phase is a more promising means of having this much mass present, but the timing 
must be right and it still would not explain the slow multiple shell systems. Slowing a massive (IOM0), fast (500 km/s) shell like 
that of r) Carinae down to 100 km/s requires 40Mq of material, and is essentially impossible. 
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individual energetic photons which we discuss below. If the dust temperature is controlled by radiative heat- 
ing, then dust can form once small grains will not be heated above the dust destruction temperature T d . If 
we consider only the mean temperature of the grains, then dust can form outside radius (IDrainell2011r) 



„ _ / L p hotQp(Tphot,a m in)\ _R* {L p hot\ ( Qp(Tphot,a m in)\ 

f ° rm ~ \l67raT d 4 Q P (T d ,a min )J ~ ~2 \~lT J \J d ) \ Q P (T d ,a min ) J ' ' 

52x W Lphot V / V 1500K V ^^ r ^^»"») V /2 cm 
Vl0 6 W V Td ) V Qe(Td,a min ) J 



(10) 



where Qp(T,a m i n ) is the Planck-averaged absorption efficiency for the smallest grains. To simplify many 
subsequent results, we define 

n _ Qp(Td,amin) nu 
Urat — n (T 7 j l 11 ,* 

SjP\lphoti a min) 

so Rform = (R*/2)(L p hot / L*) 1 ^(T* /T d ) 2 Q ra / 2 . Note that the Planck factor for the star is evaluated at the 
apparent photospheric temperature T p ] mt which may not be the same as the temperature at the stellar surface 
T*. Unless T d ~ T p h ot , the corrections for the finite size of the star are unimportant and for sufficiently small 
grains the result is independent of the g rain size because the Q oc a dependence cancels. If we use the 
graphitic models of Draine & Led (|1984h and a m j n = 0.00 l^m, the Planck average for small graphitic dusts 
is approximately a power law 

/ T \ 3/2 

Q P (T)(fim/a min ) ~ 0.42 f J (12) 

for 1000 < T < 50000 K. Transient photospheric temperatures are generally T p i mt ~ 7000 K, so the formation 
radius for graphitic dusts is approximately 

L6x10 (jfej cm - (13) 

The Planck averages for small silicate dusts cannot be reasonably approximated as a simple power law, but 
a reasonable piecewise approximation is 

log 10 [Qp(T)(fim/a min )] ~ -0.66-0.34? 3 + 1.15f 3 2 for 10 3 < T < 10 4 (14) 

where t 3 = log 10 (r/1000 K) and 

logio [Qp(T)(pm/a mi n)] ~ 0.18 + 4.90f 4 -3.21^ for 10 4 < T < 5 x 10 4 (15) 

where ?4 = log 10 (T/ 10000 K). For temperatures in the range 1000 <T d < 2000 K, the Planck averages vary 
little, so for silicate dusts and T p f, ot ~ 7000 K we find that 

R form ^ 9.8x10 (^-j cm. (16) 
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In general, including the Planck factors makes the formation radius roughly three times larger than if they 
are ignored, and a reasonably general approximation is that Rf orm ~ l0 15 (L p f !O t/l0 6 L & ) 1 / 2 cm. 

If we assume that particle nucleation occurs rapidly once the ejecta are cool enough to form dust, then 
the subsequent properties are limited by the growth of the dust particles. The collisional growth rate of a 
particle of radius a is 

da _ v c X g M 



dt I6irv w r 2 



Pbulk 



where v c is an effective collisional velocity (e.g. iKwokl I1975L iDeguchil 1 19800 . For thermal collisions, 
accreting particles of mass m a at gas temperature T, 

Ve= ^V /2 = 4. 6 r^v /2 ^V /2 «s. <i8) 



This means that growth cannot proceed by coagulation of large particles if the particle velocities are thermal 
because m a = Nmq means that v c oc N~ 1 ^ 2 and the growth rate freezes out at very tiny grain sizes. In this case, 
growth must be dominated by the accretion of monomers and very small clusters, so we can regar d m n g mn 



as eff ectively constant. Coagulation will matter if v c is controlled by turbulent motions (e.g. IVoeik et al. 



1980) with the net effect that particles can grow up to 4 times faster than by monomer accretion. The gas 
presumably cools as it expands, so we will let v c = v c o(Rf orm /R) n where Rf orm is the radius at which particle 
growth commences and n = 2/3 if the cooling is dominated by adiabatic expansion and n = 1/4 if it is 
controlled by radiative heating at constant luminosity. Other complications such as sticking probabilities 
and exhausting the condensible species can be mimicked by adjusting v c o or n. With these assumptions, we 
find that the particles grow to a maximum size of 

= VcnXgM 

arm 

* 2 5 x Iff* ( ™ Hllk \( M ( Vcn ) ( 50 ° ^ V ( 1015 Cm ^l am 
y KHM© /year J \ 0.005 J Vkm/s/ V v w J \ R form J 

VcnX g M ( U V /2 fv e \ 2 (T d \ 2 1/2 



(19) 



w i f Q ^ (20) 



~ 5 3 x lO" 3 ( MPilikQlL 2 \ / M ( 20M e \ ( v cn \ ( U \ 1/2 (vA 2 ( T d 10 4 K \ 2 
^ lO^M© /year J \ 0.005 J \ M* J Vkm/s7 \L phot ) \v w ) Vl500K J 

where v cn = v c o/(l +ri) absorbs the effects of the cooling model on a max and pbulk = Pbulk /g/cm 3 . The grain 
size grows with radius as 

R l+n 



a = a max I 1-— ^ . (21) 



Because the density is already dropping rapidly, reasonable assumptions about the temperature scaling n 
have little effect on the results. Faster cooling leads to smaller particles, but the full range from a constant 
temperature to adiabatic cooling reduces a max by less than a factor of two, and the particles are close to their 
final sizes by the time R ~ 2Rf orm . 
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If the mass loss rates are too low, then the particles cannot grow, which implies there is a minimum 



mass loss rate for dust growth of 



M > 



\6-KGM^p bldk ( 3Nm 



1/3 



£ 3.2x10 



v c x g \4npbuik 
6 /km/s\ / 0.005 



T \ 1/2 
'-'phot 



\20M & 



7* 
T d 

-•phot 



q-1/2 



1/2 



/1500K r* 



- 



T d 10 4 k 



(p 2 Nm Q ) 

Mr at 



1/3 



(22) 

year 
(23) 



where we have phrased the limit in terms of the particle number /V rather than the size a and used mo = 
mo/l2m p . Fig. Q] compares these limits from particle growth ("no growth") to those from photoionization. 
The limits from particle growth are always the more stringent. When the wind velocity has v w ~ v e , it is very 
difficult for hot stars to form dust because of the rapid increase in the wind velocity with stellar temperature, 
v w oc T*. We note that the limit at low te mperatures appears high compared to typical AGB stars (e.g. 



van Loon et al.ll2005L iMatsuura et al.ll2009f) primarily because the mass has been scaled to M* = 2OM and 
M ocM*. 

In the interstellar medium, the temperatures of the smallest dust grains are stochastic because the ab- 
sorption of individual photons can temporarily heat the grains to temperatures far higher than the equilib - 
rium temperature predicted by the ambient radiation density (e.g. iDraine & Anderson! 1 1985L iDw ek 1986). 
This effect also plays a key role in the formation of dust by transients, but seems not to have been gener- 
ally considered outside esti mates of dust formation in the colliding wind environments of WC stars (e.g. 
Cherchneff & Tielenslll995r) . Under the assumption that the ejecta must recombine in order to form dust, we 
are interested in photons with energies Eq < E < E\ ~ 13.6 eV since the hard UV photons are absorbed near 
the base of the wind. A small grain absorbing a soft UV photon will be heated well above the equilibrium 
temperature and then lose particles before r adiatively cooling. Gra ins cannot grow if this photoevaporation 
rate is faster than the collisi onal growth rate (IDraine & Salpete tl ll979h . We can estimate Eq using the models 
of iGuhathakurta & Draind (119891) for stochastic dust heating as the energy at which the grain cooling time 
scale equals the time to lose an atom from the grain. This energy depends crucially on what we view as 
the smallest number of particles in a grain because the peak temperature increases for smaller particle sizes 
and the probability of losing an atom rises exponentially with the peak temperature. If we consider a single 
Mg2Si04 molecule with N = 7 atoms as the smallest grain, then we find E\ ~ 6 eV (9.0 eV) if the grain starts 
from an equilibrium temperature of 1500 K (1000 K). If, however, we view the smallest grain as consisting 
of two such silicate units with /V = 14, then E\ ~ 13.6 eV (21.0 eV). In either case, if the transient produces 
too many soft UV photons, small grains will be destroyed by the radiation faster than they can grow. If we 
define the absorption efficiency by Q = Q'(a/X), the rate at which such photons are absorbed is 



^f^nG = Q',E /kT phot ,Ei/kT phot ) 



(24) 



where the function F is the same as for the estimate of the number of ionizing photons in Eqn. [7] but with 
G=Q' rather than G = l/x. Dust can only grow once the evaporation rate is lower than the collisional growth 
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rate (Eqn.[T7T). leading to a photoevaporation limit on the mass loss rate for dust formation of 

L phot am 



M > 



Xghc 



3mF v ^^l±£t?L f 12!M V4 (J*L-\ ( M * \ V2 (^\ (K N \ V3 ^0 (25) 
v e v c 1O 6 L V U J VlO 4 KA2OM y \X g J \ p J year' 

The factor (m^N/p) 1 / 3 ~ 2 for N ~ 7. The enormous difference between the photon and particle densities 
means that the possibility of dust formation is entirely controlled by the spectral energy distribution of the 
transient and the value of Eq. As with the recombination limits, the differences between black bodies and 
actual photospher es are crucial - the limits on M for black bodies are several orders of magnitude higher 



than those for the Castelli & Kuruca (12003b models. As we see in Fig. [T] the photoevaporation limit on M 
is a wall blocking dust formation in the quiescent state independent of mass loss rate. Thus, dust can only 
form around hot stars if they do not appear to be hot when observed from the dust formation radius. These 
stars appear to have two means of achieving this - the S Doradus phase and (great) eruptions. 

In the S Doradus phase, the luminosity of the star is little changed, L p f, ot ~ , but the stars have cooler 
photospheric temperatures, T p i wt ~ 7000 K. For our model of the S Doradus phase, we adopt the more 
common view that the lower temperature is due to a true expansion o f the stella r photosphere rather than a 



"pseudo-photosphere" formed in a dense wind (see the discussion in Vink 2009). In the S Doradus state, 
the stars have fairly high mass loss rates, M ~ 1O~ 5 -1O~ 4 M , and fast winds with ~ v e , but they cannot 
be forming significant amounts of dust even though they satisfy the condition on photospheric temperature. 
With a dust optical depth of 

Mk v ( M \ ( k v \ / 10 15 cm \ v, 

Tv 47rv w R form ~ VlO- 4 M /yeary' V 100 cm 2 /g J V R in ) v w 

the stars would become bright, hot mid-IR sources and some would be heavily enshrouded by their own 
dust, yet neither phenomenon seems to be reported. 

As we show in the middle panels of Fig.[H where we simply set the apparent photospheric temperature 
to Tphot = min(r*,7000 K), the cooler temperature is not sufficient to allow dust formation given the typical 
mass loss rates. First, the particle growth rates are too low. Second, the photosphere is still producing 
enough soft UV photons that the smallest grains still tend to photo-evaporate faster than they can grow. This 
second limit is very sensitive to the minimum photon energy Eq needed to photo-evaporate a grain, but for 
£0 = 10 eV the limit is close to the limit for any particle growth, at roughly M J> lO^MQ/year, If we lower 
the minimum energy to Eq = 7.5 eV, the required minimum mass loss rate jumps enormously because we are 
counting photons on the rapidly falling blue side of the spectrum. Thus, while the precise limits are sensitive 
to the exact choices of T p / mt and Eq, the combination of slow growth and photoevaporation mean that dust 
cannot form in the S Doradus phase. 

The final case we consider is a (giant) eruption where there is an increase in the bolometric luminosity 
Lphot > L*, the apparent temperature is cooler, T p i 10t < T*, as in the S Doradus phases, and the mass loss rates 
are much higher, M > 10~ 2 M Q /year. For sufficiently dense winds, the dust formation region sees a pseudo 



photosphere created by the non-dust opacity of the wind rather than the hot stellar photosphere. iDavidson 
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(119871) explains this as a consequence of combining a dense wind with an opacity law that is falling rapidly 
with temperature in this temperature range. Consider the Rosseland mean optical depth 



«3 



tr(R)= / P K R (p,T)dR 



(27) 



looking inwards from the radius Rj where the gas temperature is 1000 K and dust formation may be possible 
at some interior radius R. If we combine the rapidly rising p oc r~ 2 density profile of the wind with a temper- 
ature regime where the opacity rises rapidly, then there will be a tendency to produce a pseudo-photosphere 
where tr(R) = 1 near that temperature. We com puted the temperature T w ( R(tr = 1)) at the radius where tr = 1 
using the solar composition opacity models of iHelling & Lucasl (|2009l ). our standard wind density profile 
and assuming a temperature profile T w = T*{R / R*)~ l I 2 . This is not a self-consistent wind model, but the 
results are insensitive to the assumptions because the opacity and optical depth increase so rapidly towards 
smaller radii in the wind. Fig.Q]shows the consequences of using this "pseudo-photospheric" temperature in 
determining the photo-evaporation limit rather than T*, as well as the contour where T w (R(tr = 1)) = 7000 K. 

The limits now have two branches. For small M or low stellar temperatures, the wind is optically thin, 
the observed temperature is simply the photospheric temperature and the photo-evaporation limits are un- 
changed. For high M and high temperatures, the wind becomes optically thick and the observed temperature 
is of order 7000 K with re l atively weak dependencies on M and v w because of the ste ep slope of t he opacity, 
as predicted by iDavidsonl (|1987T) . As expected from the arguments summarized by IVinkl (|2009l ). the mass 
loss rates needed to form a pseudo-photosphere are higher than are typically found for the S Doradus phase. 
However, once M J> 10~ 2 - 5 M Q /year, the wind forms a pseudo-photosphere whose temperature slowly drops 
with increasing mass loss rate, which makes the photevaporation limits less sensitive to Eq than in our S 
Doradus model. Note that in both the S Doradus and eruption models the photosphere must stay in its cool 
state long enough for the ejecta to reach the dust formation radius (~ 1 year, Eqn.fTOl if dust is to form. 



Once dust forms, the size distribution then controls the opacity, 

3X d (Q\{a max )) 



^Pbulk 



where the dimensionless function 



(2a (flmax)) 



^na X !o max Q\{a)a 2 f a da 



\ amax a 3 — da 
Jo " da UU 



(28) 



(29) 



depends on the grain size distribution dn/da, the dimensionless (absorption or scattering) cross section 
Q\(a), and the fraction of the gas mass in condensed dust < X g . The function (Q\ (a max )) is proportional 
to a max for ver y small grains, where Q\ oc a and becomes constant for very large grains where Q\ becomes 
constant (e.g. iDraine & Led 1 19841) . Thus, the ratio (Q\(a max )) ja max appearing in the opacity becomes 
constant for very small gra ins, decays a s a" 1 for very large grains and has a maximum at an intermediate 



size ap ea k. At V band for a Ma this et all (| 19771 ) size distribution dn/da oc a with a range of a max /a m i n = 



50, we find a peak ~ 0.16/xm (0.45/mi) with (Q\ (a peak )) ~ 2.4 (~ 1.5) and (Q\ (a peak )) /a peak = \5Apm 1 



-12- 



(3.4/im ) for graphitic (silicate) dust. These estimates were made for the effective absorption optical depth 
(jabs(j a b s + T scat )) 1 / 2 and lead to maximum visual opacities of 

Kv,max ^ 26 °o~o7)5 cm2 /g graphitic 

~ 30 oock cm2/g silicate " (30) 

For T a bs/T scat the coefficients are 210/13 and 130/75 for graphitic/silicate dusts. In general, however, the 
size dependence of the visual opacity is relatively weak. If a max J> 1 /xm the opacity begins to drop as a~^ ax 
(Eqn.l2~TI). but this requires M ;>M Q /year, which no star seems to significantly exceed. For very small grains, 
{Qx(cimwc)) I Umax —> 7.7^m _1 (0.6^m _1 ) for graphitic (silicate) dust, so the opacity is only a factor of two 
(six) lower than the maximum opacity. If the grains cannot grow to moderate size, then the dust opacity will 
be significantly reduced. 



3. Discussion 



To summarize, the growth of dust particles in the ejecta of massive stars is limited by particle growth 
rates and photo-evaporation by soft, non-ionizing UV photons from the star. The particle growth rate is the 
limiting factor for lower temperature stars (J* <> 7000 K), and photo-evaporation is the limiting factor for 
higher temperature stars (T* <; 7000 K). While they are hot stars in their quiescent state and cannot form 
dust in their winds, the LBVs have cooler apparent temperatures, T p hot — 7000 K, in their S Doradus phases 
and during (great) eruptions. While they are cooler in the S Doradus phase, their mass loss rates are also 
not high enough for the grains to collisionally grow or to overcome the photoevaporation of small grains 
by the remaining soft UV photons. In this phase, the mass loss rates are not high enough for the non-dust 
opacity of the wind to self-shield the dust formation region, so the cooler temperat ure o f the p hotosphere 
must be due to a true expansion of the photosphere, as argued in the review by IVink! (120091) . Only in 
(great) eruptions with mass loss rates M J> lO~ 2 ' 5 M /year do these stars have the proper conditions for 
forming dust. Moreover, when the mass loss rates are this high, the wind does form a "pseudo-photosphere" 
with a temperature T p j wt ~ 7000 K that shields the dust formation region from the soft UV emission of 
the t rue stellar photosphere. This i s the characteristic "eruption" temperature of luminous blue variables 
(e.g. Humphreys & Davidsonlll994r) and it is a consequence of the steep rise of the non-dust opacity with 
temperature in this regime (|Davidsonlll987t) . 



This distinction between an expanded photosphere for the S Doradus state and a "pseudo-photosphere" 
for the eruptions is also supported by the momentum transfers needed for a radiatively accelerated wind. 
Momentum conservation requires Mvro — jL/c, where r is the n on-dust opacity source responsible for 
accelerating the wind (see iKudritzki & Pulsll2000L IPuls et al.ll2008h . For an asymptotic velocity of v^, the 
optical depth must be 

M \ / 1O 6 L \ 
M /yeary \ U ) v esc ' 



10 4 



(31) 



In the S Doradus phase, having M <, lO^MQ/year and v c 



means that r < 1 and the wind cannot form 
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Fig. 2. — Mass loss rates needed to grow dust particles to radius a max (Eqn. |2TT >. In computing the ratio of 
Planck factors Q rat we have either assumed the existence of a cooler photosphere T p j wt = min (J* , 7000 K) 
(solid) or used T p j wt = (dashed). The triangles show the estimated mass loss rates during (great) eruptions 
based on either the observed duration (filled triangles, rj Car and P Cyg) or durations estimated from the 
shell widths (open triangles). We argue in the text that these latter estimates are gross underestimates of the 
mass loss rates in eruption. The filled squares show the present day properties of the systems in Table 1. 
Objects in Table 1 noted as having exceptionally large grain sizes are circled. The temperatures are left fixed 
at the present day temperature estimates - in reality they were cooler during the eruption but we lack direct 
measurements. 
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Fig. 3. — Asymptotic expansion velocities of Type IIP SNe (top), SN "impostors" (middle) and Galactic 
LBV shells (bottom). The impostor velocities have been corrected for expansion out of the stellar potential 
well following Eqn. [33] in order to properly compare them to the LBV shells. The general pattern of the 
results is not sensitive to the details of this correction. The SNe and the shells require no corrections because 

of their high velocities and ages, respectively. The curves in the impostor and Galactic LBV panels show 

-1/3 

the expected distributions for the bins assuming a velocity-dependent rate r{v s ] ie u) oc v , for 10 km/s < 
v shell < 1500 km/s for the impostors and an observable lifetime oc l/v s heii for the Galactic shells normalized 
to the numbers of objects in each panel and excluding SN 1961V. 
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a pseudo-photosphere. In the (great) eruptions with M ;> lO~ 2 M /year, the optical depth must be r S> 1 and 
the wind must have a pseudo-photosphere in order to be radiatively accelerated. We should also note that 
once dust dust forms, it can be a significant source of acceleration because the Eddington factor for radiation 
pressure on the dust, 



_ KrpLphat ^ . ( K rp \ f L phot \ ( ' 2QM & \ 

~ 4ncGM* U V 100 cm 2 /g J \ 1O 6 L J \ M* J ' 



is large. Dust formation, as a large additional sou rce of continuum o pacity, may help to address some of the 



problems in accelerating these heavy winds (e.g. lOwocki et al.ll200 



opa 

m 



The first important consequence of this close relationship between dust formation and the need for very 
high mass loss rates is that the dust shells around luminous blue stars are formed exclusively in great erup- 
tions and so trace the history of these eruptions. Given a census of such dusty shells, their radii, expansion 
velocities, (dust) masses and optical depths in the Milky Way or other galaxies, it should be possible to 
reconstruct this dominant mass loss mechanism for these massive stars. A particularly interesting diagnostic 
is the maximum grain size. Where the total dust mass or optical depth of a shell probes the total mass lost in 
the eruption, the maximum grain size probes the mass loss rate because, as shown in Fig. [2l the maximum 
grain size is proportional to the mass loss rate, a max oc M (Eqn. [2Tb . Models of four of the Galactic shells 
appear to require a max J> 1/im (see Table 1) which strongly suggests M ;> lO~ 2 M /year or possibly even 
higher. 

A second consequence is that the mass loss rates associated with most of the shells in Table 1 are grossly 
underestimated - they are too low to make any dust let alone super-sized grains. These low estimates of M 
come from the assumption that the duration of the transient can be estimated from the radial thickness of the 
shell: At ~ AR/v s h e ii — 10 4 years since AR ~ R s heii — 1 pc and v. v / !e // ~ 70 km/s. This leads to an estimate 
of M = 10~ 4 M Q /year for M^u = M & that is not very different from many of the present day winds which 
are not making dust, as illustrated in Fig. 12 The flaw here is that the observed spread in radius probably 
comes from temporal and azimuthal variations in velocity rather than the duration of the transient, just as 
we see in rj Carinae. This is proved by the simple geometric observation that all shells have comparable 
thickness ratios, which is the characteristic of a spread in velocity: AR = Avt and R = vt so AR/R = Av/v is 
independent of time. If it were due to the duration of the transient then AR = vAt and R = vt so AR/R = At /t 
and the shells only appear geometrically thin as they become old. Roughly speaking, for every shell with 
a 2: 1 thickness ratio there should be one which is a filled sphere just finishing its eruption, and this is not 
observed. 

It should be possible to determine the geometric structure of these shells in so me detail because many 



of the central stars are known to be sig nificantly variable (e.g. 77 Carinae, see, e.g.. iFernandez-Laius et al 



a 

2C 



2009) for a full light cur ve, or. e.g.. iMartin et all 120061 for spatially resolved data; AG Car, iGroh et al 



2009; IRAS 18576+3341. iClark et al]|2009r) . You can determine both the structure of the shell and obtain 
a geometric distance to the source by mapping the time delay between the variability of the star and the 
echoes of the variability across the shell, by essent ially the same procedure as is used in reverberation 



mapp ing of quas ars (see the review by Peterson|[l993 ) or at a less involved level in studies of SN dust echoes 



(e.g. PataJ2005 ). This would complement the proper motion measurements possible for some systems (e.g. 
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Smith et all (12004a!) for rj Carinae). The optimal wavelength is probably on the blue side of the mid-IR peak, 
at 10-20/xm to maximize the sensitivity to dust temperature variations while minimizing the direct radiation 
from the star, but scattered optical or near-IR emission is another possibility if the central star is faint enough 
to allow imaging of the shell. 

With the exception of the Great Eruption of rj Carinae (250-500 km/s), the typical expansion velocities 
of the Galactic shells are only 50-100 km/s (see Table 1, Fig. [3). As we argued earlier in §2, these expansion 
velocities are unlikely to have been significantly slowed by decelerations due to sweeping up the surrounding 
interstellar medium and so must be associated with the ejection mechanism. The relatively low velocities of 
the Galactic shells mean that comparisons of the so-called "SN impostors" to LBV eruptions require detailed 



exam ination. Fig. [3] shows the expansion velocities of a sample of normal Type IIP SNe (IPoznanski et al. 



2009), the Galactic eruptions from Table 1, and the SN "impostors" from ISmith et all (|201 lh . The latter 
have been conservatively corrected to an asymptotic expansion velocity at large radius by 

' ' 2GM * (33) 



" vt 

where we used M* = 40M© and t = 14 days. While the corrections for some of the individual objects 
are sensitiv e to the choi c e of th ese parameters, the overall results are not. In this recasting of the similar 
figure from Isrnith et all ( 2011 ). we see that almost none of the impostors have velocities similar to the 
Galactic shells. We must, however, exercise care in comparing the velocity distributions of impostors and 
Galactic shells in Fig. [3] because slowly expanding shells are detectable for longer periods of time, t s h e ii oc 
I /v shell (Eqn. |2l). If the intrinsic rate of eruptions with asymptotic velocities v s heii is r(v s heli), the number of 
observable Galactic shells is cx r(v s h e ii)/v shell, independent of any other consideration such as correlations 
between v shell and M shell or completeness. 

Fig. [3] also shows a model for the velocity distributions that is consistent with both samples. We 
assumed the intrinsic rate as a function of asymptotic velocity is a pow er law, r(v shell) oc vf^n with v m ,„ < 
v shell. < v max . We excluded SN 1961V since it was probably an SN (see iKochanek et al.ll201 lL ISmith et al. 



201 lh . but otherwise ignored other ambiguities as to the nature of th e impostor sample (e.g. the very different 
physics of SN 2008S and the NGC 300-OT, see iKochaneklboilh . The best fitting model has a ~ -1/3, 
Vmin — 10 km/s and v max ~ 1500 km/s. Given the numbers of objects, the uncertainties are large (-0.75 < 
a < 0.15 for an order of magnitude change in K-S test probabilities). If the true rate is r(y 'shell) oc v^ hel! , then 
the difference can be interpreted as a velocity-dependent completeness c(v s heii) oc v°jjj^. For example, if the 
true rate is independent of v s heii W = 0), then either the impostor sample is incomplete at low velocities or 
the Galactic sample is incomplete at high velocities. There can be additional biases created by asymmetries 
in the ejection velocities, since the early time velocities may represent the fastest expanding material while 
the late time shell emission may be dominated by the slowest moving material - however, only some 50% 
of the shells are strongly aspherical (?) and the factor of ~ 2 asymmetry in rj Carinae is not large enough to 
represent a significant bias. 

We should note that velocity is not the only parameter in which there is essentially no overlap between 
the Galactic and impostor samples. First, the eruption time scales of the only two Galactic systems where 
they are known, r\ Carinae and P Cyg, are an order of magnitude (or more) longer than those of almost all 
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impostors (years to decades versus months, see lSmith et al.ll201ll for a summary), even though they are the 
only Galactic systems with relatively high velocities. Second, the ejected masses of the impostors almost 
certainly have to be far smaller than the typical Galactic shell. Assuming the impostors are radiatively 
driven, energy conservation means that the upper bounds on their ejected masses are <, O.1M , while the 
typical Galactic shell has a mass ^ 1M . Like the velocity distribution, the mass differences can be driven 
by the lifetimes of the Galactic shells, t s i ie u oc M sheU (Eqn. |2]). Nonetheless, while the Galactic sources 
and the impostors may be produced by a single process with a broad parameter range (velocity, mass, time 
scale), the observed Galactic and extragalactic sources basically sample completely different regions of that 
parameter space. If these two populations are to be unified, then the statistics and properties of the local 
systems need to be systematically determined, and the completeness of the extragalactic surveys need to be 
improved. 

In addition to (probably) being the dominant mass loss mechanism for massive stars, the eruptions may 
have a comparable importance to SNe as a source of dust, particularly as a source of large grains due to 
their favorable conditions for particle growth. Suppose that the final pre-SN mass of the stars undergoing 
eruptions is Mo and that fraction / of the lost mass is in eruptions producing material with a dust-to-gas ratio 
of X d . For SN to dominate the dust production by massive stars, they must produce 



M d , SN ^0.09 / 



X d \ /40M G 



0.005 J \M, 



1 erupt 



0.35 



1-0.26^ + 0.26- M ° 



M erupt 



M erupt 



M, 



up 



0.35 



M e 



M, 



up 



0.35 



M Q (34) 



of dust per SN. Unlike §1, w e hav e allowed M up to be fini te. We know either from the arguments of 
Humphreys & Davidson! (| 19841 ) and ISmith & Owockil (|2006ar ) or our estimates of shell statistics in § 1 that 
/ 0.5, so eruptions dominate the dust distribution unless M d SN <; O.O2/M if we are conservative (M erupt = 
4OM , M = 1OM , M up = 100M Q ) or M d ,SN £ 0.06/M© if we are more liberal {M emp t = 2OM , M = 5M , 
M up = 30 0^/^ ). Moreover, the presence of a dense circumstellar medium may also enhance dust production 
by SNe dsmith et alJl2008ah . While SN dust produc tion rates are uncert ain, few S Ne show evidence for 
producing this amount of dust (see the discussion in lMatsuura et al.ll2009h . although IMatsuura et all (1201 ll) 
subsequently reported the detection of O.4-O.7M of dust associated with SN 1987A. 

There is also increasing evidence that such high mass loss phases are crucial to understanding SNe 
on two levels. The first problem is simply one of rates. For the nominal parameters suggested by our 
discussion of the abundance of shells in §1, the rate of eruptions must be roughly the same as the rate of 
SNe. However, the rates of the faint Ty pe Iln SNe generally believed to correspon d to such transients are 
signi ficantly lower (e.g. iLi et al.l 1201 lh . This is consistent with the arguments in iThompson et all (|2009) 
and iHoriuchi et al. ( 20 lib that there is strong evidence for incompleteness in My < -16 mag tra nsients. 
Many candidate eruptions are significantly fainter, with -10.5 < My < -15 mag (ISmith et al.ll201 lh . so the 
completeness of these surveys for eruptions is presumably still worse. Local surveys need to find these 
fainter transients in order to make a complete inventory. The second issue is that there appears to be a 
mismatch between massive star formation rates and SN rates of almost a factor of two ( Horiuchi et alJboil ). 
If the underlying estimates of massive star formation rates and SNe rates are correct, then either many of 



these fainter transients need to be SNe or there must be a significant population of failed SN (Hori uchi et al 
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20111 ). This is another facet of the need to correctly classify the impostors in order to understand their 
statistics. 

The final issue we consider is that some supernova show evidence in their evolution that they are inter- 
acting with the massive dense shel ls of material created by these eruptions. The most dramatic examples are 
the hyperluminous Type Tin SNe (ISmith & McCravll2007h. a lthough thes e seem to be rela ted to low metal - 



licity environments (IKozlowski et al 



20 1( 



Stoll et all 



201 lh - However, Fox et al.1 (1201 11) argue that many 



Type Iln SNe show evidence for CSM interactions requiring the dense shells produced by eruptions. More 
generally, the existence of a dust echo from an SNe implies an eruption within 10 3 -10 4 years whenever the 
progenitor is a hot star if our theory that dust only forms in eruptions is correct!! In fact, some hyperlumi- 
nous SNe shows evidenc e for the presence of two shells - an inne r one to boost the total luminosity and 
an outer dusty shell (e.g. ISmith et alJl2008bl IKozlowski et alJl2010h . The existence of any such correlation 
has dramatic implications for the cause of stellar eruptions. In order to produce any strong SN interaction 
phenomena, the shell of material must have been produced within time t int = 10 1 5 — 10 2 5 years of the SN. 
Suppose fraction fi nt ~ 10 -10 -1 of SNe require the CSM densities of eruptions, then the time period At 
prior to the SNe over which the ejections can be occurring is 



tint ( 40M o 
At = 0AN erupt ^- — 

fint erupt 



1.35 



10 4 



erupt 



\ /0.01\ /40M G 



300 years / 



J V fint J V 



. M erU pt 



1.35 



years. 



(35) 



The existence of dust echoes leads to a similar conclusion - while tj nt is larger, their incidence in SNe 
fi„ t is higher. W hile this point has been made before in a qualitative sense (see, e.g. ISmith et al.ll201Cj 
Smith et al.ll201 lh . Eqn. [35] makes it quantitatively clear how strong a constraint results from the existence 
of any such correlation. Moreover, the parameters chosen for the scaling in Eqn. [35] may be significantly 
overestimating At. If any such correlation exists, massive shell ejections are forced to be associated with 
the very last phases of massive star evolution, roughly to the onset of carbon burning, and this suggests that 
th e underlying driving me chanism is post-carbon ignition nuclear burning instabilities (see the discussion 



in 



Smith & McCray 



2007|). Unfortunately, the only known mechanism of this kind, the pair instability SN 



dWoosley et al.ll2007h . requires very high masses (M* J> 100M Q ) and should not function at the metallicities 
of any of these nearby examples. There could s till be a strong metallicity effect because of the dependence 
of line-driven stellar winds on metallicity (see IPuls et al.l l2008b . As mass loss by normal winds becomes 
less efficient, stars may be more dependent on eruptions for mass loss, although this begs the question of 
how eruptions might become more efficient at lower metallicity. This question could be addressed by the 
investigating the statistics of LBV eruptions and shells as a function of environment. 
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2 This discussion does not apply to dust echoes from the SNe of red supergiants where it is feasible to produce the dust in a slow, 
steady wind. 



Table 1. Summary of Galactic LBVs With Dusty Shells 



Object L* r* M now 

v w,now Rsheii v sheii M s heli ^sheii a max References 
L K M /yr km/s pc km/s M & M /yr fim 



7? Car (1840) 


10 6.7 


30000 


1O-3.0 


500 


0.08 


250/500 


15 


10 00 


1 D97,S09 


77 Car (1890) 










0.03 


140/300 


0.1 


1O-2.0 




Wray 17-96 


10 63 


13000 


10 -5.5 


100 


1.0 




10 


1Q -3.0 


E02 


AG Car 


10 6.2 


29000 


10 -4.8 


110 


0.80 


70 


25 


10-2.5 


10 S91,L94,V00 


G79.29+0.46 


10 6.1 


25000 


10 -6.0 


110 


1.8 


30 


14 


10-3.3 


W96,V00b,J10 


G26.47+0.02 


10 60 


17000 


10 -4.0 


= 200 


2.3 




1.9 


10-3.4 


C03 


PCyg 


10 5.9 


19000 


10 -4.5 


190 


0.07 


136 


0.1 


10 -2.0 


N01, S06 


Wra 751 


1Q 5.8 


30000 


10-" 


= 500 


0.34 


26 


1.7 


10 -3.6 


1 H91b,deW92,V00 


IRAS 18576 


1Q 5.8 


15000 


10 -4.2 


160 


0.15 


70 


10 


10 -3.2 


U01,U05,C09,B10 


W243 


1Q 5.8 


18000 


10 -5.4 












C04 


Hen 3-519 


1Q 5.7 


28000 


10 -3.9 


365 


1.1 


61 


0.66 




S94 


HR Car 


10 5.6 


14000 


10 -5.7 


145 


0.3 


30/100 


3.0 




H91a,L96,M02 


G24.73+0.69 


10 5.6 


12000 


10-5.0 


= 200 


1.6 




0.5 


10 -4.0 


C03 


HD 168625 


10 5.4 


14000 


10 -5.9 


180 


0.48 


19 


0.25 


10 -3.7 


1 N96,P02,O03 



Note. — Luminosities and temperatures are from 


Smith et al 


(2004b 


). Other sources are BIO (Buemietal. 
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), C03 (C 


lark et al. 2003). C04 (Clark & Nesrueruela 2 


(X)4 C09 


(Clark et al. 2009). D97 (Davidson & Humnhrevs 


1997 


). deW92 Jde Winter et al. 1992). E02 (Esanetal 


2002), H91a (H 


[utsemekers & van Droml Il991ah. H91b 


(Hutsemekers & van Dromlll991bl). J10 (. 


imenez-Esteban et al. 


20ld), L96 (1 


^amers et al. 1996). L94 jLeitherer et al. 


1994), M02 (Machadoetal1l2002h. N01 (Naiarro 


200 r 


. N96 (Notaetal. 


19961). O03 ( 


O'Hara et al. 


2003), P02 


(Pasauali et al 


2002), S91 


Smith 1991 


). S94 (Smith etal. 


1994), S06 


Smith & Hartisan 2006b). S09 


(Smith 


2009D, VOOa ( 


Voors etal. 2000a). VOOb ( 


Voors et al. 


2000b). U01 (I 


Jeta et a 


. 2001), U05 ( 


Umana et al. 


2005 N 


, W96 



(IWaters et all 1 1 9961) . The ejected mass estimates are generally derived from the mid-IR dust luminosities and assume 
a dust-to-gas ratio of X j = 0.01. For 7/ Carinae the dust mass estimate agrees with estimates of the ejected gas mass 
(|Smith & Ferlandl 120071) . The S06 estimate for P Cyg is a gas mass estimate. The dust content of P Cyg is uncertain, 
although there is a mid-IR excess (see S06). 
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